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RECOIL  DISTRIBUTIONS  IN  SOME  PROTON  REACTIONS 


Single  event  upset  in  electronic  memory  devices  has  caused  a  great  deal 

of  concern  in  the  last  two  or  three  years.  ^  One  particular  area  of 

concern  is  upset  produced  in  memories  on  board  spacecraft  which  operate  in 

2 

the  earth's  radiation  belts.  Recent  NHL  calculations  have  shown  that 
upset  rates  can  be  very  high  in  spacecraft  memories.  Though  protons  do  not 
deposit  energy  sufficiently  fast  to  produce  upsets  directly,  some  of  their 
secondary  products  are  sufficiently  heavy  and  energetic  to  do  so  and  occur 
with  sufficient  frequency  to  contribute  to  observed  rates.  This  paper 
calculates  the  energy  distribution  of  the  recoil  nuclei  from  the  two  most 
important  proton  induced  reactions,  Si(p,pa)Mg  and  Si(p,2p)Al.  The  dose 
density  from  the  recoil  is  easily  obtained  from  this  and  is  briefly 
discussed.  It  is  assumed  for  each  case  that  the  two  light  products  are 
emitted  essentially  simultaneously  so  that  the  recoil  velocities  for  both 
steps  of  the  two  step  decay  are  calculated  using  the  mass  of  the  residual 
nucleus. 
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An  excited  nucleus  with  velocity  vQ  decays  first  with  a  recoil 
velocity  distribution,  D  (v  ),  spherical  relative  to  vQ,  and  again 

-V- 

with  d2^v2^'  also  spherical  relative  to  v^.  In  order  to  calculate  the 
distribution  of  the  energy  deposited  by  the  recoil,  it  is  first  necessary  to 
obtain  the  combined  distribution  of  the  two  steps  in  the  frame  of  v^. 
Agsuming  the  decayti  are  by  evaporation J 


Di(V  “  ci<Ei-Bi)e 


Vi 


i  =  1,2, 


Ei^Bi 


(1) 


where  and  t^  are  reaction  constants  and  c^  is  determined  by 


J  dE  D^( E ) , 


(2) 


The  energy  parameter  is  center  of  mass  reaction  energy  so 


„  1  2 
Ei  "  2  W 


m. 
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(VV' 


M1  =  V 
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M(  2\c) 
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Carrying  out  the  normalization,  one  obtains 


Since  it  is  our  intention  to  apply  the  vector  convolution  theorem  in 
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velocity  space  it  is  first  necessary  to  change  the  given  distributions 
from  energy  densities  to  velocity  space  densities.  Prom  (2) 


'/ 


1  3  m  /  D(v)vdv 


(  4tt) 


'lm/v'1 


D(v)d  v. 


(5) 


Thus  to  convert  D  to  a  velocity  space  density  one  writes 


D'(v) 


(4irv)  3m  D(v). 


(6) 


Assuming  no  correlation  between  and  the  combined  distribution  is 

-V  ~v 

given  as  a  function  of  u  =  v^  +  v  ^  in  the  same  frame  as  v^  by  the 
convolution  theorem. 


D,(u)=»(4tt)  2m1m2  J'J  v1iv21  Di(vi)n 

1im2  V2L  WW^l 


2(v2,5(u  -v1-v2)d3v1d3v2 


=  (4t)  m. 


,  v2=  |  u-v1|  , 


=  (4ir)“2m1m2yo1D2  dv^W^^ 

=  (87r)"imlra2/DlD2  V1V2X  dvldV 


(7) 


where  u  *  cos  9  Note  that  v2  =  0  does  not  occur.  Since 


2uv1ui  +  v£ 


V 


2 


u  - 


the  entire  dependence  is  in  v  dD  .  A^so  one  has 

holding  u  and  fixed, 


v2dv2  -  -uvld^ 

-1 ,  ,  -1, 

v2  d!^  =  -(uv  )  dv  . 

As  p^  goes  from  -1  to  1,  v2  goes  from  v+  a  u+v^  to  v 
For  the  p  integral  then  one  has 


f  V2_1  °2dPl  =  “(UV1)"1  f 


v- 


°2dV2 
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2,-1  /  ,1  2  ,  -It  ®,v2-BJ/t 

lV  J  (2m2V“B2)e  2  2  2 
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,-1  Vfc2  /  ,  2 


•(uv1)  (2)'£(t2m2)  41  e  1  j  (w“-B2/t2)e  "  dw 


-1  _t  B  /t 

(uv^)  (2t2m2>  2  e 


-w 


+  (2B2/t2-l)  /  e 
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w+ 


w+ 


=  (uv 


,-J2  B2/fc2  f  -w2  1  , 

,  )  ( 2t„m„ )  e  l  we  +  —  ( 11 )  2(  2 


(2B„/t  -l)erf (w) 


The  point  u  a  0  will  not  be  a  problem  after  changing  back  to  energy  density 
below.  Since  there  is  no  remaining  angular  dependence,  D'(u)  =  D'(u)  and, 


D'  (u) 


/S  V5  Vfc2  /\  n{  , 
[iq)  8  J  aVi(,,i)p(V 


where 


F(vx) 


=  ^we  w  +  ~  (n 


(2B2/t2-l)erf  (v 


“*  "(^f)  lutVl>'  “  ""  |  (^)*  l“'vi|  '(tf)*} 


vm  =  (2B1/ta1) 2  .  (11) 

For  the  reactions  of  interest  B2  _>  B]_  an<^  32^m2  —  Bly/inl  30 

lower  limit  for  v^  includes  that  for  v,,.  Also  in  evaluating  (11), 
whenever  (B2/t2)^  >  w+,  F,  and  thus  the  integrand,  vanish  at  this 

particular  v^. 

Of  course  at  this  point  one  has  only  the  spherical  distribution  relative 

to  v  .  This  is  a  convenient  point  at  which  to  check  normalization.  One 
o 


requires 


= J" D*  (u)d^u 


4™R  J  D’ (E)udE 


so  a  simple  check  is  just  a  sum  over 
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(13) 


D(u)  =  4TTUM'1  D'  (u) 


on  a  grid  linear  in  energy. 

In  the  numerical  evaluation  it  was  convenient  to  program  and  F  as 
statement  functions.  The  error  function  (erf)  is  a  single  precision  library 
function  and  ought  to  be  fairly  fast.  Still  each  call  of  F  is  roughly 
equivalent  to  the  evaluation  of  two  one-dimensional  integrals.  Use  of 
Simpson  rule  is  made  so  one  loop  for  evaluation  and  one  for  integration  was 
written  to  allow  printing  one  complete  integrand  for  error  checking. 
Another  convenient  error  check  was  to  require  the  sign  of  F(v^)  to  be 
positive..  The  constants  are  =  2.4,  =*  4.6,  t^  “  t?  =*  2.45  in 

MeV.  Units  of  the  are  MeV  ^  so  writing  the  masses  in  MeV  gives  the 
velocities  in  units  of  light  velocity,  c.  Since  the  integral  is  in  velocity 
and  the  normalization  check  is  in  energy,  the  corresponding  loops  are  on 
different  grids  designed  to  evaluate  overlapping  regions  in  D^tv.^)  and 
D'(u).  Thus  it  was  convenient  to  begin  each  loop  with  a  unitless  energy 
scale  using  as  the  unit.  Then,  for  example,  if  a  range  in  E(u)  from  0 
to  10  (10  B  MeV)  gives  an  adequate  evaluation  of  the  appropriate 

*3 

range  for  v, /v  would  be  1  to  10  .  While  in  principle  there  is  a 

1  m 

serious  complication  of  upper  limits  required  by  numerical  evaluation,  this 
is  ignorable  in  practice  because  of  the  property  of  these  distributions  in 
having  compact  support.  Thus  any  accuracy  desired  within  computer  precision 
may  be  attained  simply  by  extending  the  upper  limits  sufficiently. 

The  program  RECOIL  (listed  in  the  Appendix)  has  been  written,  the  first 
half  performing  the  preceding  calculations,  and  the  results  from  this  part 
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are  given  in  Tables  I  and  It.  Upper  limits  for  the  unitless  energy  grids 
were  49  for  (p,p<*)  and  36  for  (p,2p),  corresponding  to  evaporation  energies 
of  117.6  and  86.4  MeV.  For  comparison  the  peak  of  the  evaporation 
distribution  is  at  B^+t^,  or  7.05  and  4.85  MeV,  respectively.  The 
overall  sums  obtained  were  0.996  and  0.989,  respectively.  It  may  be  seen 
below  that  extending  the  energy  range  decreases  this  error  still. 

Consider  now  the  transformation  of  these  distributions  to  the  frame  of 

■4* 

v  •  Assume  absorption  of  a  30  MeV  proton  in  the  initial  excitation, 
forming  a  compound  nucleus.  One  has 


M  v  =»  M  v 
p  p  co 


1  M  2  1  2 
ttMv  =»  -t-  M  V  +  E 
2  p  p  2  c  o  x 


(14) 


and  obtains  for  the  excitation  energy 


M 

(1  -  ^-)T,  T  =  M  v?" 


M 


2  p  p 


(15) 


28 

Taking  =  M(  Si+p)  one  ob>..  .ms  E^ 


1  M  2 

77  M  V 

l  CO 


28.96  MeV,  E 


1.04,  v  ~  0.008945c.  Thus  the  transformation  is  in 

o 


part  accomplished  by  the  substitution 


2  2 
u  “V 


2w  cos  6 
o 


(16 
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Table  I 


PROGRAN  RECOIL  OH  (P, PA) 


SPHERICAL  ENERGY  DISTRiaUTIOW  CPER  HEV) 
ENERGY  QIS7RI3. 


U*GuOO  O.OOOOuU  1*6444 
G.U967  0.028580  1*7411 
0.1935  0*09*224  1.8378 
0.2902  0.173525  1.9346 
0.3869  0.246142  2.0313 
0.4836  0.302075  2.1280 
0.5804  0.343079  2.2247 
0.6771  0.372533  2.3215 
0.7738  0.393493  2.4182 
0.8706  0.40834*  2.5149 
0.9673  0.418841  2.6117 
1.0640  0.426263  2.7o84 
1.1607  0.431509  2.8051 
1.2575  0.4352*1  2.9ol8 
1.3542  0.437758  2.9986 
1.4509  0.436514  3.0953 
1.5477  0.436139  3.1920 


0.429517 

3.2888 

0.080470 

0.418171 

3.3855 

0. 069456 

0.402229 

3.4822 

8.059784 

0.382284 

3.5789 

0.05133 1 

0.359196 

3.6757 

0.043968 

fl-333924 

3.7724 

0.037581 

0*307411 

3.6691 

0.032056 

0*280503 

3.9659 

0.027292 

0.253903 

4.0626 

0.023196 

0*228177 

4.1593 

0.019682 

0*203730 

4.2560 

0.016673 

0.180844 

4.3528 

0.014104 

0.159683 

4.4495 

0.011916 

0.140330 

4.5462 

0.010053 

0.122791 

4.6430 

0.008471 

0.107026 

4.7397 

0.007129 

0.092954 

OVERALL  GUN  0.99630^' 
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Table  II 


PROGRAM  RECOIL  ON  (P.2P) 


SPHERICAL  ENERGY  DISTRIBUTION  (PEA  REV) 


ENERGY  OISTRI8. 


0.0000 

U.OOoOUO 

1.0797 

0.288428 

2.1161 

0.0432 

0.923879 

1.1228 

0.247268 

2.1593 

0.0864 

1.042541 

1 *  1664 

0.211004 

2.2025 

0. 1296 

1.979965 

1.2092 

0. 179297 

2.2457 

0.1727 

l.£  93790 

1.2524 

0.151765 

2.2889 

0.2159 

1.099337 

1.2956 

0.128003 

2.3321 

0.259  1 

1.  101681 

1.3388 

0.107606 

2.3752 

0.3U23 

1.  102705" 

1.3820 

0, *90182 

2.4184 

0.3455 

1.103173 

1.4251 

0.475367 

2.4616 

0.3887 

1.103055 

1.4683 

0.062821 

2.5048 

0.4319 

1.399393 

1.5115 

0.  052236 

2.5480 

0.4750 

1.087444 

1.5547 

0.043337 

2.5912 

0.5182 

1.063878 

-.5979 

0.035878 

2.6344 

0.5614 

1.027549 

1.6411 

0.029644 

2.6775 

0.6046 

0.979225 

1.6843 

0.U2444U 

2.7207 

G .6478 

0.920916 

1.7274 

0.020128 

2.7639 

0.6910 

0.855260 

1.7706 

0.016544 

2.8071 

0.7342 

0. 785046 

1.8138 

0.413578 

2.8503 

0.7774 

0.712890 

1.8570 

0.011127 

2.8935 

0.62U5 

0.641U4A 

,  -1.9002 

0*00910.7 

2.9367 

0.8637 

0.571315 

1.9434 

0.007444 

2.9798 

0.9059 

0.505061 

1.9866 

0.406077 

3.0230 

0«95u 1 

0.443215 

2.0298 

0.004955 

3.0662 

0.9933 

0.386352 

2.0729 

0. 004036 

3.1094 

1.0365  0.334745 


OVERALL  SUM  0.988936 


0.0&32B5 
(f  .002670 
U .002168 
0.001760 
0.001427 
0.001156 
0 .000935 
0. 000757 
0.000612 
0. 000494 
0.000399 
0.000322 
0.000259 
0.000209 
0.000168 
0.000135 
0. 000109 
0. 000038 
0.000070 
0.0 0005 6 
U.0OO045 
0.000036 
0.000029 
0.00UO23 
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~r  —r 

redefining  0  as  the  angle  between  v  and  v.  Similarly  defining  9,/  as 

Q  1 

+  -*■ 

the  angle  between  vq  and  u,  one  notes  that  (13)  has  the  form 


D  ( u ) 


ON 

dEldni 


(17) 


where  d  =  d  u  ^  even  though  it  lacks  dependence.  The 

tranformation  is  tharefora  completed  by  the  Jacobian. 


dN 

dEdfi 


3(E 

3(E,S7) 


dN 

dE  d^ 


(18) 


3(E,£2)  u 

Calculations  of  this  transformation  to  the  lab  have  been  added  to  RECOIL, 
forming  the  second  half.  In  addition  lab  solid  angle  d&  is  averaged  to 
obtain  dN/dE.  Further  the  norm  and  average  energy  are  obtained.  The 
resulting  distributions  are  in  Tables  III  and  IV  and  Figure  1.  These  tables 
show  peaks  near  1.35  and  0.97  MeV  and  average  energies  of  2.55  and  1.55  MeV 
for  (p,pa)  and  (p,2p)  respectively. 


Now  notice  that  from  (18)  one  easily  obtains  an  energy  deposition 
distribution,  given  the  recoil  energy  loss  dE/dx, 


dN  =  dE  dN 
dxdQ  dx  dEdf) 


(20) 


If  one  assumes  constant  energy  loss  (a  crude  approximation)  and  sets  dE/dx 
=  e,  then  one  recognizes  that  the  multiplication  of  (18)  by  a  constant  is 
removed  by  renormalization.  One  then  views  (20)  as  nothing  but  a  unit 
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Table  III 


PROGRAM  RtodlL  ON  CP»PA) 


SPHERICAL  ENERGY  OISTRIdJTI'IN! 

POlUTi;  Abl  lMAX  i  9.67t8  OUCtHAXJS  w.Ot'oOoG 


NORM:  6.9998a8 


LAB  ENERGY  DUUIbUlIaN  (PER  MuY) 


ENERGY 

0  «  0 

DIgTRIo. 

o.oCoOwO 

2.7 

8.20521 1 

5.4 

8 .8  5o29d 

w  •  1 

6.1a«652 

2.8 

0. 199989 

5.5 

0. 6  5288 9 

u  «  c 

6.  1  728o9 

2.9 

8  «  A  9453  9 

5.6 

6  .6  47  9a  8 

6.3 

0.  197u3a 

3.8 

8. 188982 

5.7 

0.u44u43 

w  •  4 

w  •  2 1  a  85  7 

3.1 

a. a83269 

5.8 

6.0  48  391 

C.5 

6<  At  1444 

3.2 

0.177499 

5.9 

6.0 36955 

1*0 

0.22  7  /o9 

3.3 

0.1/1649 

6 . 0 

U.8AA732 

J  .  7 

6  »  t  A  21  i  7 

3.* 

8. 165788 

6.1 

8 .6  3c  732 

o  .8 

0. A  i4  86  A 

A. 5 

8.159883 

6.2 

0.627945 

..i 

b.  A  a  t>7  5  5 

A  .  O 

8. 15a98  6 

6.3 

8.025345 

i.  j 

b . 2386  56 

3.7 

8  •  A  4  81/91 

6.4 

6.822955 

U1 

o  *  2  389 a  3 

3.8 

8.142218 

0.5 

6  .6  2o  76  5 

1  •  2 

v. t  A  94  a  3 

3.9 

8.136371 

6  .6 

0.6 18728 

l.j 

G.a3959S 

<♦.9 

w.  13u555 

6.7 

8.8 16858 

i  •  4 

6. 239563 

A.  A 

8.  124777 

o.8 

8.8 15  186 

1.5 

w.a A9AS4 

4.2 

6 • 1 a904  a 

6.9 

6.8 1 a6  3 8 

i  >6 

o  *  <138761 

4  .  A 

o ■ 11aa5o 

7  •  u 

8*6 1a2 1 1 

i  .  1 

6.237857 

‘♦.4 

8 . 10  77b  6 

7.1 

8.0  16  94  4 

1.8 

v. 2A6575 

4.5 

o  .  102 i 19 

7.2 

0.6 L 9793 

1.9 

o«c3494c 

4.6 

0  »  v  965 8  A 

7. A 

u .008  7o  7 

2  .  w 

8.23277a 

4.  7 

8 . 8  9 1 i 3o 

7.4 

8.007792 

% 

<  •  4 

b.2A.  La9 

.4.8 

o.o 85754 

7.5 

8  .800946 

t  •  £ 

8  .22o99a 

4.9 

O.b  8o484 

7.6 

J.O06177 

2.3 

0  *  22  A4b  8 

A  .8 

0.8  75342 

7.7 

6.6o54d7 

2.4 

0.217371 

5.1 

o. 878323 

7.8 

6 .80437  3 

2.5 

6.214957 

3.2 

0 . y  6  54  7  7 

7.9 

8 .804  AA  A 

4.6 

w*2  iw  242 

A  •  3 

8 .8  60  79b 

8.3 

6.603824 

lab  norm 

:  9 • 94i 59fc -6 1 

r:  AVE 

5  2 . 555u 

11 


Table  IV 


PROGRAM  RECOIL  ON  LP.2P) 


SPHERICAL  ENERGY  DISTRIBUTION: 

POINTS:  129  EMAX!  5.5278  OU<ENAX)S  0.000000  NORM:  0,988988 


LAB  ENERGY  DISTRIBUTION  (PER 


ENERGY  DI SIR  I b. 


0.0 

0.000000 

1.7 

0.1 

0.157603 

1.8 

0.2 

0.231629 

1.9 

0.3 

0.287849 

2.0 

0.4 

0.334630 

2.1 

0.5 

0.374476 

2.2 

0.6 

0.407510 

2.3 

0.7 

0.431137 

2.4 

0.8 

0.439071 

2.5 

0.9 

0.442378 

2.6 

1.0 

0.442665 

2.7 

1.1 

0.441257 

2.8 

1.2 

0.438627 

2.9 

1.3 

0.431510 

3.0 

1.4 

0.412558 

’  3.1 

1.5 

0. 393624 

3.2 

1.6 

0.373189 

3.3 

NORM 

I  9.85881E-01 

EAV6 

ME  V) 

0.352809 

3.4 

0.063937 

0.332330 

3.5 

0.053826 

0. 312182 

3.6 

0.044920 

0.292375 

3.7 

0.037175 

0.272955 

3.8 

0.030526 

0.25  3928 

3.9 

0.024877 

0.235322 

4.0 

0.020127 

0.217077 

4.1 

0.016164 

0.199204 

4.2 

0.012897 

0.181680 

4.3 

0.010234 

0.164526 

4.4 

0.008075 

0.147774 

4.5 

0.006328 

0.  131577 

4.6 

0.004927 

0.116076 

4.7 

0.003831 

0.101432 

4.8 

0.002959 

0.087793 

4.9 

0 (002266 

0.075273 

5.0 

0.001739 

1.5531 
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change  and  accepts  the  numerical  distribution  (18)  as  giving  the 
distribution  per  £  ^  Mm,  the  distance  in  which  1  MeV  is  lost.  Instead  of 
assuming  constant  energy  loss  one  may  utilize  the  code  E-DEP-l6  to  include 
straggling  effects.  A  minor  change  will  be  required  however  since  the 
energy-deposited  report  of  E-DEP-l  does  not  include  ionization  loss,  a 
significant  part  of  dose.  In  making  this  change  it  will  be  permissible  to 
ignore  radiation  loss,  a  few  percent  effect.  Short  of  using  E-DEP-l  but 
better  than  assuming  constant  dE/dx,  one  can  find  expressions  giving  the 
variation - 

In  fact  £  =*  1  MeV/hm  is  only  about  twenty  percent  low  when  averaged  over 
the  first  few  Pm.  Thus  the  implication  of  Figure  I  is  that  significant 
energy,  deposited  in  a  spherical  region  of  radius  3  to  4pm,  should  be  added 
to  that  previously  considered. 
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APPENDIX 


SOURCE  LISTING  ASC  FAST  FORTRAN  COMPILER 


CSN 

STATEMENT  CP  OPTIONS  »  <M,X) 

0001 

PROGRAM  RECOIL 

C  SET 

UP  FOR  SICP.P  ALPHAJ.  FOR  <P,2P)  SET  B2  =  Bl,  MAMP  C  CHANGE  MR 

C  TRANSFORMATION  OF  SPHERICAL  DISTRIBUTION  TO  LAB  ADDED 

0002 

IMPLICIT  REAUM) 

0003 

DIMENSION  SG(24l),DU<82>,tU(82),0V(66)*EVC66) 

000* 

DIMENSION  OV2(22)*EV2<22).OV3(22)»EV3C223 

0005 

EQUIVALENCE  COV2,OV<233).(EV2*tV<23)>.CDV3,OV(45))t(f V3.EVI45) 

0006 

DATA  MC, MP.MA, MR/931. 5016, 938. 2796, 4. 0026 *23. 98 504/, B1.B2.T 
*/2.4,4«6*2.45/«UE/1.0/»A/Q./»H»Jl»Kl/82»241»66/»HT/27.97693/ 

0007 

EUV)=<V*V*M1-B1  >*RT 

0008 

01<V)=RT*E1(V)*EXP(-F1CV)) 

0009 

FW<V)=V*EXP(-V*V)+FC*ERF(V) 

0010 

J2= Jl-1 

0011 

OVl=fl./FLOAT<J2) 

0012 

na«ma*mc 

0013 

MR=Mfl*MC 

0014 

M0=0.5*(MR4MP)/MP 

0015 

M1=M0*MR 

0016 

M2*MR*(MR+MA)*0.5/MA 

0017 

rt=i.o/t 

0018 

VW=SQMT<N2*«T) 

0019 

F-M0*VW*EXP<82*RTi 

0020 

P43 AT  AN< 1.03 

0021 

FC=<2.0*824RT-l.0)*S0RT(P4) 

0022 

VU=S0RT(  B1  /Ml'J 

0023 

Sl=2.*OVi*VU/3. 

0024 

VT’SQRTC  B2*RT) 

0025 

PRINT  61 

0026 

61 

FORMATE '  PROGRAM  RECOIL  ON  (P,PA>'> 

0027 

DU(13=0. 

0028 

EU(1)=0, 

0029 

DO  1  1=2*11 

0030 

US5QRT<DF*(I-1))*VU 

00  31 

EU(I)=0.5*MR*U*U 

0032 

00  2  J=1.J1 

00  33 

V1  =  VU*U.O+OV1*(  J-U) 

00  34 

SG<J)=0. 

00  35 

VP*< Vl*U)*VU 

0036 

IFCVP.LE.VT)  GOTO  2 

00  37 

VM»ABS<V1-U)*VW 

0038 

IF(VM.LT.VT)  VM=VT 

0039 

SGI  J)*F«<V(0-F«CVP) 

0040 

1F(SG< J3.GT.-1E-6)  GOTO  8 

0041 

PRINT  69*SG<J),V1,U,J,I 

0042 

69 

FORMATE*  ***NEGATIVE  INTEGRA NO »V1»U»J»I S'* IP 3E 11. 2*2143 

0043 

8 

SG< J)=F*01( V1)*SG< J) 

0044 

2 

CONTINUE 

0045 

S  =  0.5*<SGCJl)-SGO>) 

0046 

DO  4  J=2 . J2.2 

0047 

4 

S=S42.*SG(J)4SG< J-l) 

STATEMENT 


CP  OPTIONS  *  C*,X> 


0063 
0054 
0055 
0056 
0057 
0058 
0059 
0060 
0  061 
0062 
006  3 
0066 
0065 
0066 
0067 
0068 
0069 
0070 
0071 
0072 
0073 
0076 
0075 
0076 
0077 
0078 


0U(I)-S*S1 

AsA+nuo) 

A*A*EU<23 

P0INTS:'.I6. 

*-  £  NAX ! *  »  F7 .6  » '  DU<ENAX):*.F9.6,  NORN-  ,F9.6> 

HT*HT**C 

HE-jMB+MT 

EP=30.*NP/NE*NH/NE 

evo)=o. 

0V(1)=0. 

00  10  X=2»*l 
E  V(K )®0 » 1 *<  A-l ) 

E2*2.0*SQKT<EV<IO*EP9 

S  =  Q. 

00  11  J=1»100 

ro*i 

HU=1.0-T2*J-l)/lOQ. 

E0*£V<tO4£P-£2*NU 

Fl*SQRTCEVtO/EO> 

00  16  I*lO» 11 
IPTEO.LT ,EU(D)  GOTO  15 
CONTINUE  , 

GOTO  16 

0=OU<  IO)*(OU(I)-OUC  10>>*<£0-EIKIO»/<EUU>-EUU03> 

S3S*0*P1 

CONTINUE 

UVCK)=S*0.01 

CONTINUE 

FORMATC///'  LAB  ENERGY  DISTRIBUTION  (PER  NEV)  t 

*  PRINT  65,TEV(K).OYCK3tEV2<K)»OV2CK)»EV3<IO»OV3ClO»H=l»Z2  3 

FORNATt  3<F10.1»P10.6)) 

S=0. 

A=0. 

DO  20  K=l»*l 

s*s+nv<n 

A  =  A*TEV<K')*O.05)*OVCK) 

A-A/S 

S*S*f V<2) 

FORNATC//'*  LAB  N0RNS'*1PE13.5»'  EAVE* ' .0PF7.4 ) 

END 
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